Magnetic Compton profiles of Fe and Ni corrected by dynamical electron correlations 
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Magnetic Compton profiles (MCPs) of Ni and Fe along [111] direction have been calculated using a 
combined Density Functional and many-body theory approach. At the level of the local spin density 
approximation the theoretical MCPs does not describe correctly the experimental results around the 
zero momentum transfer. In this work we demonstrate that inclusion of electronic correlations as 
captured by Dynamical Mean Field Theory (DMFT) improves significantly the agreement between 
the theoretical and the experimental MCPs. In particular, an energy decomposition of Ni MCPs 
gives indication of spin polarization and intrinsic nature of Ni 6 eV satellite, a genuine many-body 
feature. 



PACS numbers: 71.15.Mb, 71.20.Be, 75.30.-m 

The magnetic Compton scattering is a well-established 
technique for probing the spin-dependent momentum 
densities of magnetic solids^. Compared with other 
experimental techniques, Compton scattering offers sev- 
eral advantages. Compton scattering is an inelastic pro- 
cess, in which an energetic photon collides with a sin- 
gle electron and transfers energy to it. Since the scat- 
tering is from a single-electron and (to a good approx- 
imation) occurs at a single point in space, the process 
must be incoherent and is supplying an average over real 
space. Therefore, Compton scattering is related directly 
to the electronic ground state, whereas other spectro- 
scopic methods (e.g. photoemission spectroscopy) in- 
volve excited states. 

In addition Compton scattering allows for a rather fun- 
damental test of the theories used to describe the spin de- 
pendent momentum density, since these theoretical meth- 
ods are tailored to give predictions for the ground state 
properties. Several theoretical methods have been used 
in the past to describe the electron momentum density 
and to analyze the experimental MCPs 3 -"— . Most of the 
corresponding calculations were done within the local 
spin density approximation (LSDA) for the exchange- 
correlation potential. In general, the theoretical profiles 
obtained using LSDA show an overall agreement with 
the experimental measurements, except within the region 
p z < 1 a.u.. It was pointed out that the discrepancy be- 
tween the experimental and theoretical MCP could be at- 
tributed to missing electron-electron correlations in band 
models^. Along this line, recently, several theoretical 
methods beyond LSDA have been applied in order to 
describe the features of the MCP which cannot be ex- 
plained using the LSDA-based approach and to improve 
the agreement with exp er iment£ ' 1 1 1 1 2 . The LSDA + U 
method has been applied by Tokii et al£ to calculate the 
MCPs of Fe, showing improvement for [100] and [110] di- 
rections but still underestimates near the origin the MCP 
for [111] direction. Also, one should mention the calcula- 
tions of the Ni MCP done by KuboAi implemented within 



the GW scheme based on the FLAPW method. His cal- 
culations are in overall agreement with experiment, but 
still remarkable discrepancies are found for the Ni [110] 
and [111] MCP spectra. 

The electronic structure of fee Ni has been subject of in- 
tensive studies as a prototype of itinerant electron ferro- 
magnets, since they indicate a failure of the one-electron 
theory^ 3 . - — . The LSDA calculations for fee Ni cannot 
reproduce some features of the electronic structure of 
Ni observed experimentally. The valence band photoe- 
mission spectra of N i 17 ' 18 shows a 3d-band width that 
is about 30 % narrower than obtained from the LSDA 
calculations^. Second, the spectra show a dispersion- 
less feature at about 6 eV binding energy (the so-called 6 
eV satellite) 19 ! 20 , which again cannot be reproduced by 
the LSDA calculations. Third, the magnetic exchange- 
splitting is overestimated by LSDA calculations^ com- 
pared with the experimental data^i. On the other hand, 
an improved description of correlation effects for the 3d 
electrons via the LSDA+DMFT— gives the width of 
the occupied 3d bands of Ni properly and reproduce the 
exchange splitting and the 6 eV satellite structure in the 
valence band. 

In view of these LSDA+DMFT improvement :] 22 ' 23 upon 
the magnetic properties of Fe and Ni, the comparison of 
experimental MCP with the LSDA+DMFT theory pro- 
vides some new information besides offering a test for 
the impact of electronic correlations. In particular we 
demonstrate here that the LSDA+DMFT calculations 
improves also the agreement between theory and exper- 
iment for Fe and Ni MCPs. In the following we briefly 
discuss the theoretical approach and present the calcu- 
lated MCPs of Fe and Ni together with the experimental 
data. Finally, the LSDA and LSDA+DMFT calculated 
MCPs of Ni have been decomposed and the contribution 
of different energy windows in the valence band have been 
compared in order to extract the features of electron cor- 
relations and to show their energy-dependency. Such a 
decomposition provides the evidence of the connection 
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between the MCP contribution in the lower part of the 
valence band and the existence of 6 eV satellite, both 
features being captured only within LSDA+DMFT. 
The calculations were done using the spin-polarized rela- 
tivistic Korringa-Kohn-Rostoker (SPR-KKR) method in 
the atomic sphere approximation (ASA}2&. The com- 
putational scheme is based on the KKR Green function 
formalism, which makes use of multiple scattering the- 
ory, and was recently extended to compute MCPs^Zr— . 
The spin projected momentum density n ma (p) (where 
m s =t (i) ) is computed using the LSDA(+DMFT) 
Green's functions in momentum space as 

n m M = ~ r %G^ DA[+DMFT \p,P,E)dE. 

In order to analyze the momentum density and the cor- 
responding MCPs in different energy ranges, we use a 
decomposition of the above formula in the form: 

n msAE {p) = - f 1 ^G^+ DMF ^(p,p,E)dE, 

T J Ei 

where AE = E2 — E\ represent the width of the energy 
window. The MCP seen in each energy window AE is 
obtained by performing a double integral in the momen- 
tum plane perpendicular to the scattering momentum p z 

Ji S ° A{+DMFT \p*) = J J (MP)-ndP))dp x dp y . 

Here the electron momentum density for a given spin 
orientation is given by fifm(p). The area under the 
MCP is equal to the spin moment per Wigner-Seitz cell: 

J^ M+DMFT) (Pz)dp z = ^ A(+DMFT) . In the 
actual calculations the experimental lattice parameters 
of Fe and Ni have been used (a Fe / Ni = 0.287/0.352 
nm) . The exchange-correlation potentials parameterized 
by Vosko, Wilk and Nusair— were used for the LSDA 
calculations. For integration over the Brillouin zone 
the special points method has been used 3 ^. In addition 
to the LSDA calculations, a charge and self-energy 
self-consistent scheme for correlated systems based 
on the KKR approach with the many body effects 
described by the means of dynamical mean field theory 
(DMFT) has been applied^ 3 .. As a DMFT solver 
the relativistic version of the so-called Spin-Polarized 
T-Matrix Fluctuation Exchange approximatio n 32 ! 33 
was used. The realistic multi-orbital interaction has 
been parameterized by the average screened Coulomb 
interaction U and the Hund exchange interaction J. 
The J parameter can be calculated directly within the 
LSDA and is approximately the same for all 3d elements. 
The parameter U is strongly affected by the metallic 
screening and it is estimated for the 3d metals to lie 
between 1 and 3 eV. In the present calculations we used 
J = 0.9 eV and U = 2.3 eV. 

The MCPs of Ni [111] calculated on the basis of the 
LSDA and LSDA+DMFT, respectively, are shown in 



Fig. [2 together with the experimental data. The 
Gaussian broadening applied to the calculated MCPs 
corresponds to the experimental resolution. The experi- 
mental MCPs stemming from Dixon et al*£ have been 
normalized to the experimentally determined spin mo- 
ment (0.56 /is). After broadening, the calculated KKR 
MCP spectra have been normalized to the calculated 
spin moment (0.6 /j b by LSDA and LSDA+DMFT). 
As can be seen, in the high-momentum region (p z > 2 
a.u.) the correlation effects have a small influence on 
the magnetic spin density. In the momentum region 
< p z < 2 a.u., taking into account the electron 
correlations by LSDA+DMFT approach improves the 
agreement with the experimental spectra considerably. 
Our LSDA+DMFT calculations can reproduce the dip 
in the [111] profile at ~ 0.8 a.u., which was clearly 
underestimated by the LSDA calculations. 

The LSDA and LSDA+DMFT calculated MCPs of Fc 
[111] are shown in Fig. [^together with the experimental 
spectra of McCarthy et ali^i. The experimental MCP 
has been normalized to a spin momentum of 2.07 \xb- 
The calculated spectra have been convoluted with a 
Gaussian of 0.42 a.u., corresponding to the experimental 
resolution. After convolution, the calculated MCPs 
have been scaled at a spin momentum of 2.3 [Ib (the 
LSDA calculated MCP) and 2.19^5 (the LSDA+DMFT 
calculated MCP), respectively. Although both calcu- 
lated MCPs show agreement with the experiment in 
the high momentum region, the shoulder at ~ 0.5 a.u. 
is diminished in the LSDA+DMFT calculated MCP, 
improving the agreement with the experimental spectra 
also in the low momentum region. 

The theoretical MCP of Ni [111] has been decomposed 
into contributions from different energy windows. To 
illustrate the decomposition into energy windows, the 
LSDA calculated MCP stemming from the energy bands 
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FIG. 1: MCPs of Ni [111] calculated via the KKR method 
within the LSDA and LSDA+DMFT approach. The theoret- 
ical curves have been convoluted with a Gaussian of 0.43 a.u. 
FWHM, corresponding to the experimental resolution. The 
experimental data stem from Dixon et al.— . 
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FIG. 2: MCPs of Fe [111] calculated via the KKR method 
within the LSDA and LSDA+DMFT approach. The theoret- 
ical curves have been convoluted with a Gaussian of 0.42 a.u. 
FWHM, corresponding to the experimental resolution. The 
experimental data stem from McCarthy et al.— . 
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FIG. 3: The theoretical MCPs of Ni [111] obtained from 
LSDA calculations. The MCPs have been decomposed into 
contributions of two energy windows: [-1.0 Ry, -0.4 Ry] and 
[-0.4 Ry, E_p]. 



below the Fermi level in the range [-0.4 Ry, Ep] and 
[-1.0 Ry, -0.4 Ry] are shown in Fig. [3] As can be seen, 
the main contribution of the LSDA calculated MCP is 
stemming from the energy bands situated between -0.4 
Ry and the Fermi level. The bands at energy lower 
than -0.4 Ry have just a small positive contribution 
to the MCP within the momentum range p z < 1 a.u.. 
The corresponding decomposition has been performed 
as well for the LSDA+DMFT calculated MCP and the 
results are shown in Fig. The contribution of the 
energy bands situated between -0.4 Ry and Fermi level 
has an important negative polarization for p z < 1 a.u.. 
According to the s, p and d-electron decomposition 
(not shown here) the negative contribution stem from s 
and p electrons. The overall negative polarization of s 
and p electrons is increased in the LSDA+DMFT MCP 
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FIG. 4: The theoretical MCP of Ni [111] obtained from KKR 
LSDA+DMFT calculations. The MCPs have been decom- 
posed into contributions of two energy windows: [-1.0 Ry, 
-0.4 Ry] and [-0.4 Ry, E F ]. 



compared with the LSDA approach. 
An essential feature is the important positive contribu- 
tion of the electronic states in the energy window [-1.0 
Ry, -0.4 Ry] to the MCP, that is similar for different scat- 
tering directions. The MCP carries information about 
all spin-polarized electrons in the system and about 
their localization. A broad contribution in the MCP is 
an indication for dominating localized spin states^. As 
the LSDA+DMFT MCPs are broader than the LSDA 
ones we have a clear indication that the spin magnetic 
densities have the tendency to localize in the presence 
of the electronic correlations. As this MCP broadening 
is seen in the energy range [-1.0 Ry, -0.4 Ry] where the 
well known feature of the correlated electronic bands of 
Ni, namely the 6 eV satellite (0.44 Ry) is situated, a 
direct connection between this two correlation features 
is presumable. As was shown by earlier calculations 22 
and was confirmed by photoemission experiments^, the 
6 eV satellite is spin polarized and accordingly has to 
be connected with the MCP. The general interpretation 
of the 6 eV satellite relates this feature to an excited 
state involving two 3d holes bound on the same Ni site, 
therefore it is not accessible to any LSDA calculations. 
On contrary the LSDA+DMFT approach is able to 
capture such processes via the explicit existence within 
the interacting Hamiltonian of the four index form of 
the Coulomb matrix. Although the correlation induced 
satellite is absent in Fe, the proper description of the 
angle-resolved photoemission spectra cannot be done 
by LSDA but by LSDA+DMFT^. In the case of Fe, a 
broad contribution of the LSDA+DMFT MCP is also 
obtained in the energy range [-1.0 Ry, -0.4 Ry] and 
consequently a similar tendency of localization of the 
spin density is expected. 

In conclusion, the MCPs of Fe and Ni have been 
determined using the SPR-KKR band structure method 
within the LSDA and LSDA+DMFT approach, re- 
spectively. The influence of electron correlations on 
the MCP of Fe and Ni [111] has been discussed. For 
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high transfer momenta ((p z > 2 a.u.) no signihcant 
corrections due to dynamical electronic correlations is 
seen in the Compton profile, therefore the dynamics of 
the process can be captured equally well by a LSDA 
approach. On contrary for small momentum transfer, 
a clear evidence for the interplay between the energy 
transfered to electron and the electronic correlations 
is seen: including the local but dynamical self-energy 
leads to an improved MCP spectra. In addition the 
decomposition of the Ni [111] MCPs shows a large 
and broad contribution by DMFT+LSDA approach 
stemming from the energy window between the bottom 
of the valence band and 0.4 Ry binding energy. For the 



corresponding MCP decomposition the LSDA approach 
shows just a small and narrow contribution stemming 
from the same energy range. We consider this feature as 
a consequence of the localization tendency of the spin 
density due to electronic correlations. 
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